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Abstract. Sommerfeld paradox roughly says that mathematically Couette 
linear shear is linearly stable for all Reynolds number, but experimentally ar- 
bitrarily small perturbations can induce the transition from the linear shear 
to turbulence when the Reynolds number is large enough. The main idea of 
our resolution of this paradox is to show that there is a sequence of linearly 
unstable shears which approaches the linear shear in the kinetic energy norm 
but not in the enstrophy (vorticity) norm. These oscillatory shears are single 
Fourier modes in the Fourier series of all the shears. In experiments, such lin- 
ear instabilities will manifest themselves as transient nonlinear growth leading 
to the transition from the linear shear to turbulence no matter how small the 
intitial perturbations to the linear shear are. Under the Euler dynamics, these 
oscillatory shears are steady, and cat's eye structures bifurcate from them as 
travelling waves. The 3D shears U {y, z) in a neighborhood of these oscillatory 
shears are linearly unstable too. Under the Navier-Stokes dynamics, these os- 
cillatory shears are not steady rather drifting slowly. When these oscillatory 
shears are viewed as frozen, the corresponding Orr-Sommerfeld operator has 
unstable eigenvalues which approach the corresponding inviscid eigenvalues 
when the Reynolds number tends to infinity. All the linear instabilities men- 
tioned above offer a resolution to the Sommerfeld paradox, and an initiator 
for the transition from the linear shear to turbulence. 
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1. Introduction 

The most influential paradox in fluids is the d'Alembert paradox saying that a 
body moving through water has no drag as calculated by d'Alembert via inviscid 
theory, while experiments show that there is a substantial drag on the body. The 
paradox splitted the field of fluids into two branches: 1. Hydraulics — observing 
phenomena without mathematical explanation, 2. Theoretical Fluid Mechanics — 
mathematically predicting phenomena that could not be observed. A revolutionary 
development of the boundary layer theory by Ludwig Prandtl in 1904 resolved 
the paradox by paying attention to the substantial effect of small viscosity in the 
boundary layer. Prandtl's boundary layer theory laid the foundation of modern 
unified fluid mechanics. 

Sommerfeld paradox has the potential of being the next most infiuential para- 
dox in fluids. The paradox says that the linear shear in Couette flow is linearly 
stable for all Reynolds numbers as first calculated by Sommerfeld [49j, but experi- 
ments show that any small perturbation to the linear shear leads to the transition 
to turbulence when the Reynolds number is large enough. This paradox is the key 
for understanding turbulence inside the infinite dimensional phase space. Dynami- 
cal system studies on the Navier-Stokes flow in an infinite dimensional phase space 
is still at its developing stage. In this article, we shall conduct such a study to offer 
a resolution to the Sommerfeld paradox. 

Couette flow between two parallel horizontal plates is the simplest of all clas- 
sical fluid flows with boundary layers. It is one of the most fundamental flows for 
understanding the transition to turbulence. Another basic flow is pipe Poiseuille 
flow, which is also linearly stable for any Reynolds number as shown by numerical 
computations. But Reynolds's famous experiment in 1883 showed the transition 
to turbulence for large Reynolds number. Indeed, there are lots of similar features 
in the study of transient turbulence of plane Couette flow and pipe Poiseuille flow 
[22j . We expect that some of our studies for Couette flow in this paper could be 
useful for understanding the turbulence of pipe Poiseuille flow. 

The linear stability of Couette flow was first studied by Sommerfeld in 1908 
[49] using a single Fourier mode analysis to the linearized Navier-Stokes equa- 
tions. Sommerfeld found that all eigenvalues are non-positive for all values of the 
Reynolds number, indicating the absence of exponentially growing eigen modes, 
and concluded that the linear shear is linearly stable for all Reynolds numbers. 
This fact was rigorously proved by Romanov [44 who showed that all the eigen- 
values are less than —C/R where R is the Reynolds number and C is a positive 
constant; furthermore, the linear shear is nonlincarly asymptotically stable in 
norm of vorticity. On the other hand, experimentally no matter how small the 
initial perturbation to the linear shear is, a transition to turbulence always occurs 
when the Reynolds number is large enough. 

Now we briefiy comment on some previous attempts to explain Sommerfeld 
paradox. One popular resolution, which was first suggested by Orr [39j (see also 
[51j ). is to use the non-normality of the linearized Navier-Stokes operator to get 
algebraic growth of perturbations before their final decay. (Note: non-normality 
refers to operators with non-orthogonal eigenfunctions.) However, it is not clear how 
such linear algebraic growth relates to the nonlinear dynamics (see [55]). Moreover, 
the non-normality theory cannot explain many coherent structures observed in the 
transient turbulence. Most recently, another resolution of the Sommerfeld paradox 
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from an interesting new perspective was proposed [23j by investigating the half- 
infinity Couette channel {x > 0) and e^°"^ (cr > 0) type basic modes. However, this 
does not answer the paradox when x is in the whole line or periodic setting which 
is commonly used in the literature for numerical or experimental studies. 

We believe that norms play a fundamental role in resolving the paradox. Per- 
turbations with large norm of vorticity may still have small norm of velocity, 
i.e., small energy. In experimental or numerical studies, such perturbations of 
small energy are still considered to be small, but they are outside the vorticity's 
neighborhood of the linear shear where Romanov's nonlinear stability result is valid 
[44j . Therefore, such perturbations have the potential of being linearly unstable 
and initiating the transition to turbulence. 

The main idea of our resolution is to show the existence of a sequence of linearly 
unstable shears which approach the linear shear in the velocity variable but not in 
the vorticity variable. These shears are the single modes of the Fourier series of 
all the 2D shears of the Couette flow y + X^m^i "^"j sin{my). More precisely, our 
sequence of oscillatory shears has the form 



As n — !■ oo, the oscillatory shears approach the linear shear, i.e. Un{y) ^ y in L'^ 
and L°°. On the other hand, in the vorticity variable, the oscillatory shears do not 
approach the linear shear since dyUn{y) = 1 -f 4A7rcos(4n7rj/) 7^ 1 in any Lebesgue 
norm. Thus in the velocity variable, the oscillatory shears can be viewed as the 
linear shear plus small noises. For any large n, we prove that Un{y) is linearly un- 
stable for both inviscid and slightly viscous fluids. More precisely, in Theorems 13.21 
and 14.11 it is shown that these shears are linearly (exponentially) unstable for both 
Euler equations and Navier-Stokes equations with large Reynolds numbers. Under 
the Navier-Stokes dynamics, our shears are not steady rather drifting slowly. By 
viewing them as frozen, we investigate the spectra of the corresponding linearized 
Navier-Stokes operator (Orr-Sommerfeld operator). Moreover, numerical simula- 
tions [25j indicate that its unstable growth rate does not depend on n substantially, 
implying that as n — s- 00, its unstable growth rate does not shrink to zero rather 
approach a positive number. Such a linear instability will generate the transition to 
turbulence no matter how small the initial perturbation added to the linear shear 
is, as long as the Reynolds number is large enough to realize the linear instability. 
Numerical simulations [25] indicates that initial perturbations to Un{y) can lead to 
a transient nonlinear growth of 10 times. 

In Theorem l5.ll we prove the bifurcation to nontrivial travelling solutions to 2D 
Euler equation, near the oscillatory shear C/„(y). The streamlines of these travelling 
waves have the structure of Kelvin's cat's eyes. In Theorem 16. 1[ we also show that 
3D shears {U {y, z) , 0, 0) in a neighborhood iW^'^ {p > 2) in the velocity variable) 
of any linearly unstable 2D shear (including our Un{y)) are linearly unstable too. 

In recent years, there has been a renaissance in numerical dynamical system 
studies on fluids. The focus was upon three classical flows: plane Couette flow, plane 
Poiseuille flow, and pipe Poiseuille flow. Basic flows of them are the linear shear 
for plane Couette flow, and the parabolic shear for both plane and pipe Poiseuille 
flows. Original studies on transition to turbulence for these flows were conducted 
by the Orszag group [40j [41) . Primitive steak-roll- wave coherent structures were 
discovered. The Orszag group also emphasized the importance of the slowly drifting 
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states like our Un{y)- In 1990, Nagata made a breakthrough by discovermg 3D 
steady states (fixed points) in plane Couette flow [38 . There are two branches of 
3D steady states, called upper and lower branches. The same 3D steady states also 
appeared in plane Poiseuille flow 56^. Most recently the same 3D steady states 
were also discovered in pipe Poiseuille flow and they become traveling waves with 
constant speeds in the laboratory [9] [22] . The pipe discovery made laboratory 
experiment easy to implement. Indeed, Hof et al. observed such traveling waves 
in experiments [19] . A steady state (fixed point) study is generally the starting 
point of a full dynamical system study which is still in progress |21j |13j |14j [53] 
[47j |45j . The most intriguing new discovery on the steady states is that the lower 
branch steady states of the three classical flows share some universal feature, i.e. 
they all approach 3D shears as the Reynolds number approaches infinity '57' [54J 
[56j . Such 3D shears {U{y, z), 0,0) are neutrally stable under the linearized Euler 
dynamics. So they are not our 3D shears in the neighborhood of our Un{y). On the 
other hand, any 3D shear {U{y, z), 0, 0) is a steady state of the 3D Euler equations; 
so it is interesting to understand the special features of the limiting 3D shear of 
the lower branch. A necessary condition on such a limiting 3D shear was also 
discovered [29j . In numerical and experimental studies on transition to turbulence, 
one often observes the steak-roll-wave structure. The streak and roll are the 0-th 
Fourier mode, while the wave is the higher Fourier mode in the Fourier series of the 
velocity field. When plotting these modes separately, the steak-roll-wave structure 
should generically be observed. Taking into account the linear instability of our 
3D shears in a neighborhood of our Un{y), such steak- roll- wave structures will be 
generated too. 

The 3D steady state argument on the Sommerfeld paradox claims that the 
stable manifolds of the 3D steady states may get very close (they cannot be con- 
nected) to the linear shear in the infinite dimensional phase space, which leads to 
the transition from the linear shear to turbulence [13 . The main criticism of this 
argument is that these 3D steady states are not close to the linear shear in the phase 
space even in the kinetic energy norm; while experiments show that no matter how 
small the initial perturbation to the linear shear is, a transition always occurs when 
the Reynolds number is large enough. We believe that the linear instability of our 
Un{y) (and 3D shears in its neighborhood) offers a better explanation for the ini- 
tiation of transition from the linear shear to turbulence. In other words, the key 
point for answering the Sommerfeld paradox is whether or not there is a linear in- 
stability happening arbitrarily close to the linear shear. Our Un{y) (and 3D shears 
in its neighborhood) does the job. Again, they approach the linear shear in the 
velocity variable but not the vorticity variable, in consistency with Romanov's non- 
linear stability theorem [44]. On the other hand, stable manifolds of the 3D steady 
states can often only be established in higher Sobolev spaces in which Romanov's 
nonlinear stability theorem prohibits them to get close to the linear shear. 

Explorations on two dimensional viscous steady states turn out to be not suc- 
cessful so far |5] [lOj . That is, the counterpart of the 3D upper or lower branch 
steady state has not been found in 2D. On the other hand, numerics shows that 
transitions still occur from the linear shear to turbulence in 2D. This further con- 
firms that the stable manifolds of the 3D steady states are not the initiators for 
the transition. Under the 2D Euler dynamics. Theorem 15.11 shows the existence of 
2D inviscid steady states with cat's eye structures. The neighborhood of these eye 
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structures might be a good place for a future numerical search of 2D viscous steady 
states. 

A rigorous proof on the existence of 3D steady states in the three classical flows 
remains open. The main difficulty is the lack of a proper bifurcation point. There 
is no hope directly from the linear shear end since it is linearly stable. From the 
infinite Reynolds number limiting shear of the lower branch, it is promising and 
challenging since we do not know much of its property except a necessary condition 
[29j . However, with the discovery of more and more 3D N-S steady states (travelling 
waves) numerically besides the upper and lower branches [15] , it is unclear if any of 
these 3D steady states are fundamentally important to the development of transient 
turbulence. In Theorem 15.11 we show that nontrivial 2D steady states bifurcate 
from the unstable shears (C/„ (y) , 0) under the 2D Euler dynamics. Since C/„ (y) 
is arbitrarily close to the linear shear in the velocity variable, so we get inviscid 
steady states near the linear shear. These inviscid steady states could provide a 
natural starting point for constructing steady states of Navier-Stokes equation with 
large Reynolds number. This is a problem of bifurcation from infinity or asymptotic 
bifurcation theory as suggested by Yudovich 59j. Of course, to get N-S travelling 
waves from Euler travelling waves, the key and difficult issue is to understand the 
boundary layers. 

Proving the existence of chaos (turbulence) for Couette flow is a far more 
difficult problem than resolving the Sommcrfcld paradox. The boundary layer 
effect adds tremendously to the difficulty. The Kolmogorov flow (periodic boundary 
condition in every spatial dimension) is a much easier mathematical problem in this 
aspect. The boundary layer is not present. Some progress on the dynamical system 
studies on the Kolmogorov flow has been made [36] [3] [4j [28] [17] [26]. So far, 
proving the existence of chaos in partial differential equations is only successful for 
simpler systems, for a survey, see [27]. 

The article is organized as follows: In section 2, we will discuss some possible 
phase spaces in which dynamical system studies can be conducted on the Couette 
flow. In section 3, we prove the inviscid linear instability of our sequence of os- 
cillatory shears Un{y)- In section 4, we prove the viscous linear instability of our 
sequence of oscillatory shears Un{y) when viewed frozen. In section 5, we prove a 
bifurcation of our oscillatory shears Un{y) to Kelvin's cat's eyes under the 2D Euler 
dynamics. In section 6, we prove the inviscid linear instability of 3D shears U{y, z) 
in a neighborhood of our oscillatory shears Un{y)- 



We are interested in fluid flows between two inflnite horizontal planes (Figure 
[T]) where the upper plane moves with unit velocity and the lower plate is flxed. The 
dynamics of such a fluid flow is governed by the Navier-Stokes (NS) equations 



defined in the spatial domain Dqo = M x [0, 1] x M, where u = (1*1,^2,^3) is the 
velocity, p is the pressure, and e is the inverse of the Reynolds number e = 1/R. 
The following boundary condition identifies the specific fiow 



2. Mathematical Formulation 



(2.1) 



ut + u ■ Vu 



Vj3 + eAu , V • 7i ; 



(2.2) 



ui(t, x, 0, z) = 0, ui(i, X, 1, z) = 1, 
Ui{t, X, 0, z) = Ui{t, X, 1, z) = 0, [i = 2, 3). 
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Figure 1. Coucttc flow. 



The linear shear is given by 

(2.3) ui =y, U2= U3 = 0. 

One can choose the infinite dimensional phase space to be 



S 



e iJioc(^oo), (s > 3), V • M 0, together with 



where H^^^ is a local Sobolev space. This phase space S is too large and too difficult 
to analyze. One can first study its invariant subspace by posing extra periodic 
boundary condition along x and z directions with periods Li and L3. Denote by 
D the partially periodic domain D = [0, Li] x [0, 1] x [0, L3], the partially periodic 
invariant subspace is given by 



e H%D), (s > 3), V • M = 0, together with lU 



(2.4) S = 

Often we are interested in the two-dimensional reduction U3 — 83 — in which case 
the phase space S is simplified further. Moreover, in the two-dimensional case, the 
flow is nicer since both 2D NS and 2D Euler equations are globally well-posed. 

Remark 2.1. By the change of variables 

the new variable v — {vi,V2,V3) satisfles the Dirichlet boundary condition at y = 
0, 1; thus 

-t-00 

V = Vn{x, z) sinmry, 



and Vn{x,z) is periodic in x and z. By this representation, we see that the phase 
space S' is a Banach manifold. In fact, the Navier-Stokes equations (|2.ip can be 
re- written in terms of the v variable as 

(2.5) Vt+v-Vv^~Vp + eAv + f , V ■ v = ; 

where / = ~ydxV — V2 (1,0,0). In terms of v, the phase spaces will be Banach 
spaces. 

In our exploration inside the phase space, we will focus not only on the neigh- 
borhood of the linear shear ()2.3|) . but also on the neighborhood of the following 
sequence of oscillatory shears 

(2.6) wi = [/„(?/) = y+ ^sin(4n7ry), < A < , U2 = "a = 0, 
which will be proved later to be linearly unstable under the 2D Euler flow. 
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All the 2D shears together form an invariant submanifold 

ui = U{y), U2 = U3 = 0, U{0) = 0, C/(l) = 1 



(2.7) A = lyes' 

where U{y) is an arbitrary function. Inside the invariant submanifold A, the dy- 
namics is governed by 

dtU^ ed^U. 
The fixed point of this equation is given by 

d^U^O, i.e. U = ciy + C2. 

The boundary conditions U{0) = and U{1) — 1 imply that 

which is the linear shear (|2.3p . That is, when e 7^ 0, the linear shear is the only fixed 
point inside A. Of course, when e = (Euler flow), A is an equilibrium manifold. 

Let 

U^y + V, 

then 



utv — td^ 

Thus 



dtV = ed^V, V{0) = V{1) = 0. 



+ 00 

V a„e^'("^)'* sinnTry, 

n=l 

where a„'s are constants. Finally the orbits inside A is given by 

-t-00 

U = y + Y^ a„e-'("^)'* sinnvry. 

n=l 

3. Inviscid Linear Instability of the Sequence of Oscillatory Shears 

Let {U{y), 0) be a steady shear of the 2D Navier-Stokes or Euler equation, e.g. 
U{y) = y. Introducing the stream function tp: 

dtp dtp 
oy ox 
and linearizing at the steady shear in the form 

^ = 0(y)e*"=^+^* = (/)(y)e'"(="-'=*) , A = -lac; 
one obtains the so-called Orr-Sommerfeld equation 

(3.1) —[52 - + [/'> -{U- c)[d^ - a^]cf> = 0, 
with the boundary conditions 

(3.2) ^ = ^' = 0, aty = 0, 1. 

For Euler equation, e = and p.ip is reduced to the Rayleigh equation 

(3.3) U"(f>^{U-c)[d^-a^]<j) = 0, 
with the boundary conditions 

(3.4) (f>^0, at y = 0, 1. 

We call (a, c) with a > 0, c G C an eigenmode of Orr-Sommerfeld or Rayleigh equa- 
tions, if for such (a, c) the equation (|3.ip or p.3p is solvable and the corresponding 
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solution is called an eigenfunction. The eigenmode with Imc > is called unstable, 
with Imc = is called neutral. As shown in last section, all the shears {U{y),0) 
are steady under the 2D Euler dynamics; while only the linear shear U{y) = y is 
steady under the 2D Navier-Stokes dynamics. Nevertheless, when the Reynolds 
number R is large (e is small), the shears {U{y),0) only drift slowly. When they 
are viewed frozen (or artificial forces make them steady) , the spectra of the corre- 
sponding Orr-Sommerfeld operator are still significant in predicting their transient 
instabilities. 

First indicated by Sommcrfeld '49', the linear shear U (y) = y is linearly stable 
for all Reynolds number. This fact was rigorously proved by Romanov [44| who 
showed that the eigenvalues of the linearized Navier-Stokes Operator satisfies 

ReA < -C/R; 

and moreover, the linear shear is nonlinearly stable for all Reynolds number in 
M^^'^([0, 1] X R^). When e = 0, for the linear shear, the Rayleigh equation reduces 
to 

(U - c)ip = 0, where (f =^ [dy - a^](/), 
which has only continuous spectrum [11) 

c e [min U (y), max U{y)] = [0, 1], 
and the corresponding quasi-eigenfunction for c = Ulyo) is given by 

If = 6{y~yo). 

In the (j) variable 

[dy-a'^](p = S{y-yo), 

that is, (j) is the Green function. In fact, the linear shear is also nonlinearly stable 
under the 2D Euler flow. In the vorticity form, the 2D Euler equation is given by 

dfl 

— +u-Vn = 0. V • u = 0; 
ot 

which has the invariants J F{il)dxdy for any F. Near the linear shear, fl = 1 + uj, 
and Lo satisfies 

duj 

— + u • = 0, V • u = 0; 
ot 

where u is the velocity corresponding to Q. Then 

J F{uj)dxdy 

are also invariant for any F. Thus the linear shear is nonlinearly stable in norm 
of vorticity for any p G [1, +cxd]. 

When e = 0, each point in A is a fixed point. The linear spectrum of these 
shear fiows has been studied a lot since Lord Rayleigh in 1880s |43j . If the profile of 
the shear does not contain any inflection point, then by Rayleigh's criterion there 
is no unstable eigenvalue and the spectrum consists of only continuous spectrum 
given by the imaginary axis. However, the existence of an inflection point is only 
necessary for linear instability and the results on sufficient conditions for instability 
remain very limited [8] |31) . To construct unstable shears near the linear shear, we 
use the following instability criterion for monotone shear flows. 



SOMMERFELD PARADOX 



9 



Lemma 3.1. Consider a monotone shear profile U{y) G C (0, 1) with inflec- 
tion points. Let y'^ , ■ ■ ■ be all the inflection points and {U'' — U (2/*)}'_o be the 
inflection values. Define 

(3.5) Qdy) - ^"^^^ 



U{y) - 

and the Sturm- Liouville operator 

(3.6) L,^ = -ip" + Q.,{y)^ 

with the Dirichlet boundary condition ^p{Q) — <yj(l) = 0. (i) If for some < i < I, 
the operator Li has a negative eigenvalue, then the Rayleigh equation US. 3\) has 
unstable eigenmodes for some intervals of wave numbers. Any end point of the 
unstable intervals is such that —aj. being a negative eigenvalue of some operator Li. 
(a) If Li > for alio < i < I, then the shear profile U (y) is linearly stable for any 
wave number. 

Lemma l3. II was rigorously proved in [33) . The following key observation is due 
to ToUmien [52' in 1930 (see also and [30]): if —a^ is a negative eigenvalue of 
Li and (ps is the eigenfunction, then (c, a) = (U^, as) is a neutral eigenmode to the 
Rayleigh equation (|3.3p . with 0^ being the eigenfunction; then one can try to find 
unstable modes (c, a, 0) (Im c > 0) to (j3.3p near such neutral mode. A variational 
formula for dc/da can be derived at the neutral mode. 

Notice that any point in A has the representation 

+ 00 

(3.7) U (y) = y + ^ sin niry, 

n=l 

SO the sequence of oscillatory shears (jl.ip can be viewed as single modes of the 
above expansion. The following theorem shows the inviscid linear instability of the 
oscillatory shears (jl.ip . Since |C/„ (y) — 2/|^oo(o i) ^ this shows that there exist 
unstable shears in any small (velocity) neighborhood of the linear shear (|2.3p . In 
particular, these unstable shears are arbitrarily small kinetic energy perturbations 
to the linear shear. 

Theorem 3.2. Under the 2D Euler dynamics, the oscillatory shears C/„ defined 
by ULl]) are linearly unstable. More precisely, there exists an unstable eigenmode 
curve (a,c{a)) with Im c(q!) > of the Rayleigh equation ^3.S^) with Un, stemming 
from a neutral mode (ia„,l/2j where an > cqu (co > is independent ofn). The 
corresponding unstable eigenf unctions are in C°° (0, 1). 

Proof. First notice that the oscillatory shears (jl.ip are monotone, since > 
0. So to show the linear instability of t/„, by Lemma [3.11 it suffices to prove that 
for at least one inflection point, the Sturm-Liouville operator (|3.6p has a negative 
eigenvalue. We choose the inflection point of J7„ at y ^ 1/2. Define 

U"{y) _ -16TT^nAsm{4mTy) 
^ ^ ^ [/(?/) -C/(i)-y-i + ^sin(4n7r2;) 
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and the Sturm-Liouville operator L — — + Q{y)- By Rayleigh-Riesz principle, 
the smallest eigenvalue Ai of L is given by 



(3.9) 



Ai = min 



dy 



Thus, to prove that Ai is negative, all we need to do is to show that the fraction on 
the right hand side of above is negative for some specific test function ip £ Hq {0,1). 
Let y' = y-\,y' & [-\, 5] and denote Q{y') = Q{y' + i). When y' € 
we have 

2 

— |4n7ry'| < | sin(4n7ry')| < |4n7ri/'|. 



— —1 



TT 



thus 



~ -leTT^nA |sin(4n7r2/')l -Wtt'^uA^ \Am:y'\ 

Q{y ) ^ ~, ~ ~ ^ 



^ |sin(4n7rj;') 



^ |4n^y'| 



-2(87r)^»M 
1 + AttA 



We choose the test function 93 as follows (Figure [2]) 

v{y') 



j_ -J_<i,'<J_- 

4n ' Sn — » — 8n' 



8n y ' 8ri ^ W — 8n , 
8n ^ w ' 8n — y ^ Sn ' 

0, -2>\y'\>i,- 

Since both Q{y') and (f{y') are even in y', we only need to estimate over the interval 





l/4n 






-3 -2 -1 
8n Sn 8n 


1 2 3 

8n 8n 8n 



Figure 2. The test function ^p. 



[0, ^] . First we have 
Notice that for r/ e [0, 



since J7„ is monotonically increasing, and 



' sm 



4n7r ( \- 1] 

An 



Ann 



1 

An 
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Thus 



Notice also that 



so 



4n 



v) > Q{y' - ^ + ^) > 0. 



On the other hand, 



dr] < 0. 



1 + 4Att 
AAtt^ 1 
1 + AAtt An 



1.] 1. 

An J 8n 

7r(l-(5) 1 
' 2-S 4^' 



where 4 An =1 — (5,0<(5<l/2for the osciUatory shears We also have 



1 

An 



1 

8^ 



1 

An 



Gathering all the above estimates together, we get 



dy' 



/o ^'^dy 
By p.9p . this means 

(3.10) 



< 



Ai < 



-liAnf 



2-5 



< 



7r(l -(5) 



- 1 



when < S < ^ K 0.533. Since < (5 < 1/2 for the oscillatory shears C/„ (y) f^ . 
the corresponding Ai is negative. Thus by Lemma 13. H we get unstable modes 
(a, c(a)) (Imc(a) > 0) to the Rayleigh equation (13. 3|) associated with J7„ (y) , for 
certain intervals of wave number a. In particular, the unstable wave number inter- 
vals include a neighborhood of a„ = %/— Xj" > CQn and for those wave numbers a, 
the unstable eigenvalue c (a) is close to the inflection value ^. Here, Cg is a positive 
constant independent of n (see p. 101) ). Since J7„ (y) S C°° (0, 1), the regularity of 
the unstable eigenfunction follows easily from the Rayleigh equation (|3.3p . This 
completes the proof of the theorem. □ 



Remark 3.3. 1). Even though they approach the Couette flow (|2.3p in the 
energy norm and L°° norm of velocity as n — > cx), the oscillatory shears (j2.6[) do 
not approach the linear shear in iJ" (s > 1) as n ^ oo; by simply noticing that 

?7,'j 1 + ATTAcos{Amry). 

Thus in our phase space S, the oscillatory shears are of a finite distance away from 
the linear shear for large n. On the other hand, in fluid experiments, it is the L°° 
norm of velocity that is observed; thus for large n, the oscillatory shears (|2.6p can 
be considered to be small noises by fluid experimentalists. 
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2) . The form of the unstable shears J7„ (y) ()l.ip is not unique. For example, 
we could choose 

(3.11) Un{v)^y + -W{nv) 

n 

where W (y) G (0, 1) is a 1— periodic function of the similar shape as siny. By 
choosing constant A properly, the proof of Theorem 13.21 can still go through to 
get the linear instability of [/„ (y). Nevertheless, we consider the oscillatory shears 
p.6p as good representatives of linearly unstable shears near the linear shear since 
they are single Fourier modes of (|3.7p . 

We also notice that the oscillatory structure of [/„ (II. ip is in some sense nec- 
essary for instability. First, for any given function / (y) G C^(0, 1), the shear 
flow Ue (y) = y + e/ {y) is spectrally stable, i.e., there is no unstable solution to 
the Rayleigh equation (13. 3|) . when e is small enough. This can be seen as follows. 
When e is small, Ue (y) is monotone. So by Lemma [3.11 (ii), Ue (y) is spectrally 
stable if all the operators Li — — ^ + Qi{y) are nonnegative, where 

'^'^'^^'uAy)-UAy.) 

and i = 1, • • • are all the inflection points of / (y). When e is small, IQiiy)]^^ 
is small so £i > which proves the linear stability of Ue (y) ■ Second, for shears of 
the form (|3.1ip with W (y) G (R) , to preserve the Couette boundary conditions 
p.2p for all large n, the function W (y) must have infinitely many zeros. So it is 
natural to choose as a periodic function with zeros. 

3) . For the Sturm-Liouville operator L — — ^ + Q{y) with Q{y) defined 
by (|3.8p . the second eigenvalue A2 > (see Appendix for a proof). Thus Ai is 
the only negative eigenvalue of L, so by Sturm-Liouville theory Ai is simple and 
the corresponding eigenfunction </)„ (y) can be chosen such that <^„ (y) > when 
ye (0, 1). This fact will be used later in the proof of Theorem 15. II 



Let D = [0, Li] X [0, 1], where Li = ^ and a is an unstable wave number for 
Un (y) with n large. As a corollary of Theorem 13.2) we have 

Theorem 3.4. The spectra of the linear Euler operator at the oscillatory shear 
1^ m H^iD) (s>3) are as follows: 

(1) There are J > 1 unstable eigenvalues and J > 1 stable eigenvalues. 

(2) The imaginary axis is the absolutely continuous spectrum. 



The first claim is proved in Theorem 13. 21 and the proof of the second claim can 
be found in 

For any linearly unstable shear flow of Euler equation, the nonlinear instability 
can be established too (i.e. [l], [16j . [32j). Moreover, the existence of unstable 
and stable manifolds near unstable shear flows was recently proved [34] for Euler 
equation. 

4. Viscous Linear Instability of the Sequence of Oscillatory Shears 

A more precise statement of the Sommerfeld paradox is as follows: 

• Mathematically, the linear shear is linearly and nonlinear stable for all 
Reynolds number R, in fact, all the eigenvalues of the Orr-Sommerfeld 
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operator satisfy the bound A < —C/R where C is a positive constant 

m- 

• Experimentally, for any R > 360 (where i? = ^ in our setting '2^), there 
exists a threshold amplitude of perturbations, of order 0{R^^) where 
1 < M < X depends on the type of the perturbations [24] , which leads to 
transition to turbulence. 

A mathematically more precise re-statement of this experimental claim is as 
follows: For any fixed amplitude of perturbations to the linear shear, when R 
is sufficiently large, transition to turbulence occurs. For any fixed R, when the 
amplitude of perturbations is sufficiently large, transition to turbulence occurs. 
There may even be an asymptotic relation between such amplitude threshold and 
R. 

Our main idea of the resolution is as follows: The oscillatory shears (j2.6p are 
perturbations of the linear shear. As n — > oo, they approach the linear shear in 

norm of velocity. They are linearly unstable under the 2D Euler dynamics. 
As shown later on, this will lead to the existence of an unstable solution of the 
Orr-Sommerfeld equation (|3.ip with U given by the oscillatory shears (|2.6p . when 
the Reynolds number R is sufficiently large. Notice that these oscillatory shears 
are not fixed points anymore under the Navier-Stokes dynamics. Nevertheless, they 
only drift very slowly. The important fact is that here the unstable eigenvalue of 
the Orr-Sommerfeld operator is order 0(1) with respect to e = as e ^ 0+. 
This fact should lead to a transient nonlinear growth near the oscillatory shears 
(and the linear shear) which manifests as a transition to turbulence. This has been 
confirmed numerically [25] . Here the amplitude of the perturbation from the linear 
shear will be measured by the deviation of the oscillatory shears from the linear 
shear and the perturbation on top of the oscillatory shears. One final note is that 
here the turbulence is often transient, i.e. with a finite life time after which the 
fiow re-laminates back to the linear shear. 

As mentioned in the Introduction, due to the recent discovery of 3D steady 
states, there has been a conjecture on the explanation of the Sommerfeld paradox 
using the stable manifolds of the 3D steady states. The stable manifolds can not be 
connected to the linear shear since it is linearly stable. The conjecture is that the 
stable manifolds can get close to the linear shear, which leads to the transition from 
the linear shear to turbulence |13) . The main criticism on this argument is that 
these 3D steady states are quite far away from the linear shear even in the velocity 
variable. For each of such 3D steady states, its stable manifold will have a fixed 
distance from the linear shear due to its Romanov's nonlinear stability theorem 
144' . When the Reynolds number R is large enough, this distance will be much 
bigger than the threshold ©(i?^'') of perturbations given above. This contradiction 
shows that the stable manifold explanation is not so satisfactory. 

Based upon the fact that the oscillatory shears (|2.6|) are linearly unstable under 
the 2D Euler dynamics, it is natural to search for an unstable solution of the Orr- 
Sommerfeld equation (|3.ip with U given by (jl.ip . The main difficulty in this search 
naturally lies at the boundary conditions. For Rayleigh equation (|3.3p . only two 
boundary conditions are necessary while four are required for Orr-Sommerfeld. 
When e is small, of course the other two boundary conditions will create boundary 
layer effects. By using the asymptotic expansion theory developed by W. Wasow 
[58] 60 years ago, we can prove the following. 
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Theorem 4.1. Let (a^,c^) he on the unstable eigenmode curve of the Rayleigh 
equation associated with the oscillatory shears C/„ il-l]) . which stems from the neu- 
tral mode (aQ, 1/2). Then for e sufficiently small, there exists an unstable eigenmode 
(a'^, c* ) with Imc* > of the Orr-Sommerfeld equation 113.1]) with U given by il.l]) . 
Whene^Q+, c* c° . 

Proof. Rewrite the Orr-Sommerfeld equation 

(4.1) [d'y-a' 

where Un is given by the oscillatory shears (|2.6p , 

j^=iaR, 6(2/) = -([/„- c). 

Let (a^jC^) be on the unstable eigenvalue curve of the linear 2D Eulcr operator 
(e = in (|3.ip ). stemming from the neutral mode (a„, 1/2) as identified by Theorem 
c" > (c" is near 1/2). We will study the segment where c^ 

a circular disc of radius 



' + [b{y) {d'y 



3.ip in the following form: 

- K] = 0, 



is sufficiently close to 1/2, and the region c G Bs^{c'^) 
< Si < and centered at c". See Figure [3l 




Figure 3. A c-plane illustration. 



We know that /7„(y) — 1/2 has only one zero y = 1/2 on the real segment y G [0, 1]. 
Thus we can find a rod region S around [0, 1], in which Un{y) — 1/2 has only one 
zero dX y — 1/2. See Figure |4l Then by Rouche theorem, when c is sufficiently 
close to 1/2, Un{y) — c also has the only zero y near 1/2: 



1 



1 + 4:TtA 



1 

c 

2 



O 



1 

c 

2 



and 
(4.2) 



K{y) = {l + 4TTA) + 0(c 



^0. 
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Define 



Qiy) 



y^—b{s)ds. 



As y circles around y once in counter-clockwise direction, a,ig{'yQ{y)) increases 
by Stt. So there are three curves Cj (j = 1,2,3) stemming from y, on which 



Re{'yQ{y)) = 0. By the derivative formula (|4.2p . there is a disc D centered at 
y — 1/2, the radius of which is independent of c, such that inside the disc, the three 
curves Cj can be well approximated by using 



Un 



{l + A7TA){y-y). 



We can choose the rod region S so thin that it penetrates D as shown in Figure 
[S] Direct calculation reveals that Ci has an angle of 7r/6, and the angles between 
neighboring Cj's are 2tt/3. When c = 1/2, along the two pieces of the interval 
[0,1] outside the disc D, arg(7Q(j/)) is not changing. Thus when c is sufficiently 
close to 1/2 but not equal to 1/2, the entire real interval [0, 1] lies inside the two 
sectors 5*1 U 5*2 {y not included. Figure [4]) . Inside the two sectors 5*1 U 5*2 , Wasow 
obtained the following results [58j: There are four linearly independent solutions 
to the Orr-Sommerfeld equation (|4.ip in the forms. 



(4.3) = e^^'-y^ 



N 



,r!=0 



, (j = l,2); 



(4.4) (^e^My) + f{Si'JS2h-^ (^ = 3,4); 
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Figure 5. Another y-plane illustration. 



where Re{'jQ{y)) < in Sj for and changes sign when y crosses into the other 
sector; Tpjn are analytic in S*! U £'2, iV is a large number, f{Si U S2) denotes any 
function which, together with all its derivatives in ?/, is uniformly bounded in 7 in 
every closed subdomain of S'l US'2, and ipiiu) s-re two linearly independent solutions 
to the e = Orr-Sommerfeld equation (|3.ip . In particular, here 4>j {j = 1,2,3,4) 
are four linearly independent analytic solutions on the entire real interval [0,1]. 
The eigenvalues of the Orr-Sommerfeld operator (|4.ip are given by the zeros of the 
determinant 

MO) 02(0) 

01(0) 0^(0) 

0l(l) 02(1) 

0i(l) 0^(1) 



A(c,7) = 



03(0) 

03(0) 
03(1) 
0^,(1) 



04(0) 

0^(0) 

04(1) 

0U1) 



Each entry of A(c, 7) is analytic in c e Bs-^{c'^), (c" 7^ 1/2). Wasow obtained the 
following expression 



A(c, 7) = 72e'^(Q(0)+Q(i)) [KAoic) + 0{j-^)] 

where 

A (r) - ^3^0) ^4(0) 
^3(1) 04(1) 

and Re{jQ{0)) > 0, Re{'yQ{l)) > 0, X is a non-zero constant. Let 

Ai(c,7)=ifAo(c)+0(7-^). 
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Notice that the unstable eigenvalue c° of the linear 2D Euler operator is a zero 
of Ao(c). Here we fix a = a^. By a proper choice of Si, Ao(c) is non-zero on 
the boundary of -65^(0"). Then when 7 is sufficiently large, by Rouche theorem. 



Ai(c, 7) also has a zero c* near c*^. As 7 — > 00, c* 



The proof is complete. □ 



Remark 4.2. The key point in the above proof is to show how to embed the 
entire real interval [0, 1] inside the interior of the two sectors 6*1 U 5*2. This key 
point was missing in |37| . 

Next we will develop an expression for the eigenfunction needed for later studies. 
First we need a lemma. 

Lemma 4.3. Re{'yQ{y)) is a strictly monotone function on y ^ [0, 1]. 

Proof. 



dy 



Re{iQ{y)) = ^[RehQiy)) ~ Re^Qim 



dy 



rv ^ 

'oiR-rR^ I + i{Un 
dy Jo 

^ aRRe\/ Ci + i{Un — c^). 



Cr)ds 



For the above to be zero, we need Ci < and Un 
the above is never zero This proves the lemma. 



Cr — 0- Since Ci > in our case, 

□ 



The eigenfunction of the Orr-Sommerfeld operator corresponding to the unsta- 
ble eigenvalue given by the above theorem is given by 

01(2;) 02(2/) hiy) (l>A{y) 
0i(o) 0^(0) ^'M 04(0) 

0l(l) 02(1) 03(1) 04(1) 

01(1) 0^(1) 0^(1) 0U1) 

where 0j (j = 1, 2, 3, 4) are given by (|4.3p - (|4.4p . We need to know the asymptotic 
property of as 7 —> 00. For this, the only trouble maker is the exponent jQ{y)- 
We choose the convention 



(4.5) 



i?e(7Q(0)) < < i?e(7Q(l)). 



(4.6) 



01 (y) 


02 (y) 





0^(0) 


0i(l) 





01(1) 





03(1) 


04(1) 


0^,(1) 


0U1) 


0i(l) 


04(1) 


0'i(l) 


0U1) 



Some of the entries in the expression of are exponentially small in 7. Dropping 
these entries, we have 

03 (y) 04 (y) 
03(0) 0UO) 

03(1) 04(1) 

0^(1) 0U1) 

0^,(0) 0^0) 

0l(l) 03(1) 04(1) 

0U1) 0^(1) 0U1) 

0l(l) 03(1) 
0'l(l) 03(1) 

The dominant terms are the third and the fourth terms, which are of order 



02 (0)01 (y) 
02(o)03(y) - 



'2(y) 
02(O)04(y) 



O (72e^'=('yQ(i))-«'=(')'Q(o))') 
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under the convention (|4.5|) . By rescaling 0, 

=7-2e«-(7Q(O))-fle(7Q(l))0. 

we see that 

0. = C[04(l)03(2/) - 03(l)04(y)] + 0(7~')- 

Notice that the quantity (j)4{l)(f>3{y) — (/)3(l)(/)4(?;) is not zero, but smaU of order 
0(7^^). It depends upon c*. When 7 ^ 00, it does approach zero. Taking 
derivatives of (j) in y, the balance of orders shifts to the first two terms in (|4.6|) . One 
can clearly see that 

U4h^ ^0{\jr'), as 7 ^00. 
In fact, it is a matter of scaling. If we scale 4> as 

for any N ^ then 

~0(|7|^-^+i), as 7 -00. 



5. A Bifurcation to Kelvin's Cat's Eyes 

The bifurcation we are seeking will stem out of the oscillatory shear (12. 6|) . 
Since there is no external parameter in 2D Euler equation, the bifurcation here is 
with regard to an internal parameter: the wave number along x-direction. The 
proof of Theorem 13.21 shows that there is a neutral mode (a, c, (/)) = (Q;„,^,0n) 
to the Rayleigh equation (|3.3p associated with the oscillatory shear [/„ {y) (jl.ip . 
Here, — = Ai is the negative eigenvalue of the operator L = —-§^7 + Q{y) with 
Q{y) defined by (|3.8p . and 0„ is the corresponding eigenfunction. By the proof 
of Theorem 13.21 > cqti. By Remark 13.31 3). we can choose (f>n{y) such that 
(pn (y) > when ye (0, 1). The following theorem shows that we can get nontrivial 
inviscid travelling waves bifurcating from the above neutral modes. 

Theorem 5.1. There is a local bifurcation curve of travelling wave solutions 
to the 2D Euler equation with the stream function = 'ip°'{a [x — ,?/), which 
stems out of the oscillatory shear t/„ {y) i2.6\) . Here a is near an > cqti (Theorem 

V) e C^'^ ((0, 27r) X (0, 1)) (0 < /3 < 1) 

is periodic and even in ^ of period 2tt, and ip" is constant on {y = 0} and {y — 1}. 
Near y — 1/2, the streamlines of these new fixed points have a Kelvin cat's eye 
structure, with a leading order expression given by \5.5\) . 

Proof. Let ^ = ax where a will serve as the bifurcation parameter. Let 
V'rcilCi y) to be the relative stream function in the reference frame {x — ^t, y), that 
is, 

i/'roi(C,y) = i'ii^y) - -jV- 

The travelling waves solves the Poisson's equation 

(5.1) a) = a -—- + -—- - /(VAci) = 
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for some function /, with the boundary conditions that tAioI takes constant values on 
{y — 0} and {y = 1}. Since we seek nontrivial traveUing waves near the osciUatory 
shear (|2.6p . we demand that the osciUatory shear (|2.6p satisfies (15.11) too. That is 

(5-2) ^ - m.,) 

where 

1 / A 



V'rci (y) = o U - o - TTIT cos (47i7r?/) 



is the relative stream function for the the oscillatory shear (|2.6p . It follows from 
(EH) that 



(5.3) mci)- ,,™i 

where (5(2/) is defined before in (|3.8p . Since both ■0*^^ (y) and (5(?/) are symmetric 
to the hne y ~ ^, we only need to check (|5.3p for < j/ < i. The function 7/)*^,) (y) 
is monotone on [O, 5] since V'*oi (?/) = ^^"(y) — 5 < on [0, 5). Therefore we have 

which determines the relation / in (15. ip . We define 

1 / 1^ ^ 

•/"(f-y) = Aci{£.,y) ^ 2 

and reduce (|5.ip to solve the equation 

with the homogeneous boundary conditions 

0(^,0) = <^(^,1) = 0. 
For /3 e (0, 1) , define the spaces 



0) = 1) =0, 27r — periodic and even in ^ 
and 

L) = y) e C°'''([0, 27r] X [0, 1]), 27r - periodic and even in ^} . 
Consider the mapping 

F((/),a2) : BxR+ ^ D 

defined by 

Then the travelling wave solutions satisfy the equation F{(j),a'^) = 0. The trivial 
solutions corresponding to the oscillatory shears (|2.6p have 

A A 
{y) = --^ cos (4n7ry) 



4n27r 471^77 
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Let — be the negative eigenvalue of — |^ + Q{y) and (pniu) the corresponding 
positive eigenfunction. The hnearized operator of F around (0*, a^) has the form 

C F^,{<j)^,al) = (^Ig^ ^dy^ /'(V-rci) 

Then by Remark 13.31 3) , the kernel of £ : _B 1— > _D is given by 

keiiC) = {(f>n{y) cos£_} , 

In particular, the dimension of ker(£) is 1. Since C is self-adjoint, (j)n{y)cos^ ^ 
R{C) - the range of C In fact, again by Remark l3.3l 3). 

Aun{B/R{C)} = 1. 

Notice that da2d^F{(t),a'^) is continuous and 

9„29^i^(0,,a^) (0„(y)cosO = ^ [0n(y) cos^] = -0„(y) cos^ ^R{L). 

Therefore by the Crandall-Rabinowitz local bifurcation theorem [6] , there is a local 
bifurcating curve (/?)) of a^) = 0, which intersects the trivial curve 



a^ ) at = at 



(f>iP) = (f>*iy) + (y) cos e + o(/5) , 



and a^(/3) is a continuous function, a^(0) = 0;^- So the relative stream function 
has the form 

(5.5) ^liy-l) +My) + PMy)cosC + oip). 

Since 4>n{y) > in (0,1), so near the inflection point y = 1/2 of the oscillatory 
shear (|l.ip . the streamlines of these travelling waves have a cat's eye structure (see 
[8]) with saddle points near (i,27rj) (j € Z). The proof is complete. □ 

Remark 5.2. 1) The small travelling waves constructed in Theorem 15. ll has a 
a;— period of the order O (i) . In particular, the cat's eyes near the oscillatory shear 
Un (y) has the spatial scale ^. 

2) Since the oscillatory shears C/„ (y) (|2.6p are arbitrarily close to the linear 
shear in norm of velocity, the travelling waves constructed in Theorem l5.1l appear 
in an arbitrarily small (L^— velocity) neighborhood of the linear shear. We note that 
there might not exist nontrivial travelling waves near the linear shear in a stronger 
norm in lieu of Romanov's nonlinear stability theorem [44j . This conjecture is 
partly supported by the following rough argument. Any travelling wave of 2D 
Euler equation satisfies the Poisson's equation 

(5.6) -Aij = g{^) 

for some function 5 in = (0, L) x (0,1), where L is the x-period and ip is the 
relative stream function. Taking x derivative of (|5.6p . we get 

-Aipx = g' (ip) ipx- 
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Note that ipx = on the boundaries {y — 0} and {y — 1}. Muhiplying above by 
Tpx and integration by parts in f2, we get 

/ / \'^il)x\^ dxdy ^ / / g' (%l>)\%l)x\^ dxdy. 
J Jn J Jn 

If the travelling wave is close to the linear shear in a strong norm (e.g. C^-vorticity), 
g should be close to 1 in norm and thus \g' (V')Iloo is very small. We have 



|V?Aa;| dxdy <\g' {tp)\^^ I I dxdy 



n 



But we also have 



iVipxl^ dxdy ^ {-Aipx^ipx) > J J lijjxl^ dxdy. 

Since = on dfl and the operator —A with Dirichlet boundary conditions on 
80, has the lowest eigenvalue tt^. So if \g' (V')Il=o < n^, we must have ipx = in 
and the travelling wave is a trivial shear flow. 

6. Inviscid Linear Instability of 3D Shears 

In this section, we show that the instability of oscillatory shears J7„ (y) (jl.ip 
persists under the 3D setting. Consider a 3D shear flow uq = {U {y, z) , 0, 0), which 
is a steady solution of 3D Euler equation for any profile U {y, z). The fluid domain 
is 

5^3 — {(x, y,z) \ Q < y < \, Lx and periodic in x and z} , 

where Lx and Lz are to be determined later. The linearized 3D Euler equations 
near mq are 

(6.1) dtu + Uux + vUy + wUz = —Px, 

(6.2) dtv + Uvx = -Py, dtw + Uwx = -P., 

(6.3) Ux + Vy + Wz ^ 0, 
with the boundary conditions 

(6.4) v{x,0,z) = v{x,l,z)^0. 

Here (m, u,w) and P are perturbations of the velocity and pressure. Consider a 
normal mode solution ^'"^('■^ (^u, v, w) {y, z) to the linearized equation, with a — 
^IT — 1'2,---). By eliminating u and P, we end up with a PDE system 

(6.5) {U - C) [Vyy - a'^V + Wyz) - UyyV - UyzW - UzWy + UyWz = 0, 

(6.6) {U - C) (Wzz - + Vyz) - UzzW - UyzV - UyVz + U zV y = 0, 

with the boundary conditions t; (0, z) = u (1, z) = 0. When U depends only on y 
and w = Q, the system (j6.5p - (16.61) is reduced to the Rayleigh equation (|3.3p for 
2D shears. So far, the understanding of the instability of such 3D shears remains 
very limited due to the complicated nature of (|6.5p - (|6.6p . Our next theorem shows 
instability of 3D shears close to an unstable 2D shear. 
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Theorem 6.1. Let Uo (y) e (0, 1) be such that the Rayleigh equation \3.3\) 
has an unstable solution with (ao,co) (ao,Imco > 0). Fixed Lz > 0, consider 
U{y,z) g ((0,1) X [QjLz)), Ly-pcriodic in z and 

(6.7) C/(l,z) = [/o(l),C/(0,z)-[/o(0). 

// ||?7 (y, z) — f/o (2/)|lvKi.p((o i)x(o L^)) (-P > 2) is small enough, then there exists 
an unstable solution e*"°^^~'^*-' (u, w, w, P) (y, z) to the linearized equation around 
{U {y, z) ,0,0) with \c—co\ small. Moreover, ifU{y,z) G C°°, then {u,v,w, P) £ 
C°°. 

The proof of Theorem 16.11 is divided into several steps. First, we give a new 
formulation of linearized growing modes for 3D shears. Fixed a = ag. Consider a 
growing mode solution ^^'^oi^-ct) j-^^ (Imc > 0) to the linearized Euler 

equation around {U {y, z) , 0, 0). Instead of studying (|6.5|) - (|6.6[) . we reformulate the 
problem in the following way. From (I6.ip - (l6.4p . we have 



(6.8) iao {U — c)u + vUy + wUz — —iaoP, 

(6.9) iao {U — c)v — —Py, ia^ {U — c)w ~ —Pz, 

(6.10) iuQU + Vy + Wz =0, 
with the boundary conditions 

(6.11) u(0,z) = i;(l,z) = 0. 
Denote — {{y, z) \ < y < 1, Lz periodic in z}, 

8^2 = {y = 0} U {y = 1} , 

V2 = {dy,dz) , A2 = dyy + dzz, 

and U2 = {v, w). We claim that: 

(6.12) / / " '^y^^ = 0, 

and 

(6.13) * U2-dl^ f U2-dl = 0. 

J{y=0} ■J{V=1} 

The identity (|6.12p follows by integrating (|6.10p in ^2 with the boundary condition 
(|g?7|) . The identity (pTTO)) is a resuh of ([Q]) . (p?7)) and the assumption that Imc > 0. 
Define 

(6.14) uj — Wy — Vz- 
Then by equation (|6.9p . 

(6.15) {U -c)uj + UyW-UzV^O. 
Taking {y, z) divergence of (|6.9|) and using (j6.10|) . we get 

(6.16) - A2P = al(U - c)u + iao {UyV + Uzw) . 
From (|6.1ip and (|6.9p . P satisfies the Neumann boundary conditions 

P,(0,z) = P,(l,z) = 
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on dfl2- Denote 

(^^2) = {feL' (n-,) \ I I f dydz - 



H^'^ (n^) ^ e H^n^) I J J^^ fdydz = 0\ 

and Q : — + to be the projector to the mean zero space L^'°. For any 
/ € L^'", denote h — (— A2)^^ / to be the unique solution in ($12) of the 
Neumann problem 

-A2h = /, in Q2 
dh 

— — = on dfl2- 
on 

Thus from (|6.16p . we get 

(6.17) QP = (-As)^' [al {U-c)u + iao [UyV + U,w)] . 

Denote B = (-A2)w^ Q = Q{-^2)n^ Q, then B : {VL2) -> H^'° {^2) is a self- 
adjoint, bounded and nonnegative operator. We rewrite (j6.17p as 

(6.18) QP = alB iUu) - ca^Bu + iaoB {UyV + U^w) . 
Multiplying (|6.8p by Q and using the equation (|6.18p . we get 



u 



(6.19) -iaoc{l + alB) 
= -lao (alB + Q) (Uu) + {a^B - Q) {UyV + U,w ) , 

where we use the property Qu = u due to (|6.12p . Let A = —iaoc, then Re A > 0. 
From (|6.19p and (|6.15p . we get 

(6.20) A ( ^ ) = ^^aoU ( ^ 
(1 + alBy^ [tao (1 - Q) (Uu) + {a^B - Q) {UyV + U,w )] 

— IUq {UyW — UzV) 

—ianU I " 
V ^ 

/ /^^ Uu dydz + (1 + alB)-' {a^B - Q) {UyV + U.t 
~iao [UyW - Uzv) 



since 



(l + alB) '(l-Q)(C/ii) 



1 r r . . \ 2.^-1 _ 1 



-r Uu dydz (1 + a^B) 1 = ^ / / Uu dydz. 

Lz J JSI2 / Lz J Jri-2 

Thus the growing mode problem is reduced to study the unstable spectrum of the 
operator A ^ T + IC, where T is the —iaoU multiplying operator and 



^ / ' I ' ' -iao {UyW - Uzv) 
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In the above definition, (Wjv) G (iJ^ ($12))^ is uniquely determined from {u,lu) G 
(L^ (^^2))^ by solving equations 

(6.21) Vy + Wz — —iaoQu, 

with the zero circulation condition (|6.13p and the zero normal velocity condition 
V = on 9ri2-This is guaranteed by Lemma |7. II in Appendix. 
We study properties of A in the next lemma. 

Lemma 6.2. (i) A : (i^ (^^2))^ {L^ (^2))^ is compact perturbation of T. 
(a) The essential spectrum of A is i [a^ min t/, ap min U] . 

Proof, (ii) is a corollary of (i) because of Weyl's theorem (see [17] or [20j ) and 
the fact that the operator J- is bounded, skew-adjoint and has the essential spec- 
trum i [ao min [/, ao min U]. To show (i), we need to prove that /C : (L^ (^^2)) 
(i^ (^^2))^ is compact. By Lemma [7.11 we have 

I ] = ||u2 11^1^2) - '^o(||divu2 |Il2 + ||curlu2 11^2) 

< Co(ao||w IIl2 + 11 ujWl^). 
Thus the linear mapping {u,uj) — > {v,w) is compact in (L^ ■ Since 

{l + alBy\B, Q 

are bounded, this proves the compactness of /C. □ 

As a corollary of the above lemma, any eigenvalue A of ^ with Re A > is a 
discrete eigenvalue with finite multiplicity. 

By using the formulation (|6.20p . the proof of persistence of instability for 3D 
shears is similar to that in |12j for 2D Euler equation. Given A with Re A > and 
U (y, z) e ((10, 1) X (0, Lz)), we define the operator 

M (A,f/) := {X-Ty^K. 

m (L2 (O2)) . Then the growing mode equation is reduced to solve M 
u \ 

for some A with Re A > 0. 

u: J 

Lemma 6.3. Consider X G C+ = {Re A > 0} and U {y,z) satisfying conditions 
in Theorem \ 6.1\ Then (i) A4{X,U) : (L^ (5I2)) (i^ (172)) is compact and 
analytical in X E C+. (ii) M. (A, C/) depends continuously on U in the following 
sense. Let V (y, z) satisfy the same conditions ofU {y, z) as in Theorem \6.1\ . Then 
for any b > 0, there exists another constant C > such that 

(6.22) sup \\M (A, U)-M (A, V)\\^(,^.,^.... < C \\U - . 

Proof. Since J- is skew-adjoint, for any A with Re A > 0, (A — J^) ^ is bounded 
and M is compact in (L^ (^^2)) by Lemma |6^ Since !F generates an unitary group 
and 

/.+00 
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SO Ai (A, U) is analytic in A in the half-plane C+. To show (ii), we write 
M{\U)- M{\V) 

= {\-T (A, U)y^ K. (A, U)-{\-T (A, V)y^ K. (A, V) 

(A, U)y^ {T (A, U)-T (A, T/)) {\~ T (A, F))"^ /C (A, U) 

+ {\-T{\V))-\lC{\U)-K{\V)) 
= 1 + 11. 



When ReA > 6, 



< 



j^. Both IC{X,U) and J^{X,U) are norm 

continuous to t/ in the norm ||C/||vi/i,pi which we show below. We use C for a 
generic constant. First, 

\\J^{X,U)-J^{X,V)\\ <ao\\U-V\\^^ <C\\U-V\\ 
Second, for any {u,uj) G (L^ (^2))^ ,we have 



1^1. p 



(/C(A,[/)-/C(A,y)) 

<c(^\\U-V\\^.^\\u\\ 
<c(^\U-V\\^^\\u\\ 
<C\\U-V\\^,^, 



L2 



L2 



L2 

|v(c/-v)IL. 
|v(c/-y)|ii. 



V 

w 

V 

w 



_2p_ 
Lp-2 



L2 



and thus 



|!(/C(A,f/)-/C(A,F))|l^. <C|1C/-1/|1^^,„. 
So both / and // above are uniformly controlled by \\U — VWw^^p- 



This proves 
□ 



The proof of Theorem 16.11 uses the following Lemma of Steinberg [48j . 

Lemma 6.4. Let T{X; s) be a family 0/ compact operators on a Banach space 
analytic in X and jointly continuous in (A; s) for each (A; s) e Ax S , where A is an 
open set in C and S is an interval in R. // for each s there exists a X such that 
I — r(A; s) is invertible, then (/ — ^(A; s)) ^ is meromorphic in X for each s and 
the poles of {I — T(A; s)) depend continuously on s and can appear or disappear 
only at the boundary of A or at infinity. 

Proof of Theorem 16.11 By definition, A is an unstable eigenvalue of A if 
and only if 1 is an eigenvalue of Ai (A, U). For Re A > and < s < 1, define 

r (A, s) = {l-s)M (A, Uo) + sM (A, U) . 

By Lemma l6.31 these operators are compact, analytic in A and continuous in s. The 
assumption on Uq (y) implies that Aq = —iaoCo is a pole of [I — T (A, 0)) ^ in the 
right half plane. Since the poles of (/ — T (A, 0)) are discrete, we can choose so 
small such that on the circle F = {A | |A — Ao| < £o}j the inverse {I — T (A, 0))~ 
exists. By the continuity property ()6.22p . 

Iir (A, s)-T (A, 0)11 < 5 \\M (A, Uo) - M {X,U)\\ < C \\U - Uo\\y^^,, . 
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Thus when ||C/ — ?7o||vt'i,p is sufficiently smah, {I — T (A, s)) also exists on the cir- 
cle r for all s E [0, 1]. Thus by Lemma [6^ there exists a pole Ai of (/ — T (A, 1)) ^ = 
{I — M (A, U))" within the disk {|A — Ao| < Eq}. For Eq small, Ai also has positive 
real part and is an unstable eigenvalue of the operator A associated with U {y, z). 
Let (m, uj) e (^L? (O2)) be the corresponding eigenfunction and (w, w) G (i/^ (^^2)) 
is determined by {u,lo) as in the definition of tC. Define 

c = — — P ~ — : — [ictQ {U — c)u + vUy + wUz) ■ 
—iao —iceo 

We check that {u,v,w, P) solves the growing mode equations (|6.8p - (|6.1ip . First, 
an integration of the equation Au = Xiu implies that J J^^ u dydz — 0, namely 
Qu = u. So (|6.10p follows from (|6.2ip . The equations (|6.8p and (|6.1ip are already 
implied in our construction. To check (j6.9p . first we note that the vector field 
iao (U — c) (v , w) is curl free because of the equation (|6.15p which follows from 
Aid — \\Lo. So there exists a scalar function P' {y, z) such that 

iao (U-c) {v,w) = -{P'^,P'y) . 

Taking divergence of above and using (I6.10p . we can set 

P' - (-Aa)^' {al {U-c)u + iao [UyV + U,w)) 

by modulating a constant. It remains to show (|6.17p . from which (|6.9p follows. We 
note that Au = Aiu implies (|6.19p . Combining with Qu — u, we get (j6.18p . an 
equivalent form of (|6.17p . 

To show the regularity of the growing mode, we look at the eigenfunction 
equation 

The eigenfunction {u,ijj) e (L^ (^2))^ implies that {v,w) G [H^ (^2))^ and thus 
/C ^ ^ ^ e (i?^ (^2))^ ■ Since (Ai — J-)^^ is regularity preserving, we have G 

(i/i(f72)) . If U e C°°, we can repeat this process to deduce that {u,uj) G 

{H''{n2)f for any fc > 1 and therefore {u,u;) G (C°° (f^s))^ Then u,v,w,P G 
C°° (fJ2). This finishes the proof of Theorem O □ 

Remark 6.5. In Theorem 16. 1[ the smallness of ||C/ (y, z) — t/o (2/)I1h'1'P(02) 
{p > 2) is required to show persistence of instability. For 2D shears, i.e. U = U (y), 
we can show that: if \\U (y) — Uo (y)llHi(o 1) i^ small enough, then linear instability 
of Uo (y) implies that of U (y). This is an improvement over the result in [12j . 
where the smallness in norm is required to prove persistence of instability for 
2D Euler fiows without hyperbolic points. The proof is very similar to the 3D case, 
although much simpler. So we only sketch it briefiy. In the 2D case, we reformulate 
the growing mode problem in terms of u. Using the notations as before, from the 
linearized equations we derive 

Xu = -iaoUu + iao / Uu dy + {l + alB)~^ {alB - Q) (Uyv) . 
Jo 
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Here, Q is the projector from (0, 1) to mean zero space and B — Q (^—-^^ Q 

where ^— ^jp-^ is the mean zero solution operator of the ID Neumann problem. The 
mapping u —i- v is defined by solving the ODE 

dv 

— = -iaoQu, v{0)^v (1) = 0. 
ay 

The rest of the proof is the same as in the 3D case, so we skip it. Smallness of 
\\U (y) — Uq (?/)||j^i(q is required in the proof because of the Sobolev embedding 
(0, 1) ^ L°° (0, 1) . 

7. Appendix 

Proof of the claim A2 > in Remark 13.31 3): By the Sturm-Liouville the- 
orem [60 , the linear operator L, 

Lip = -ip" + Q{y)^, ,^(0) = (^(1) = 

has a sequence of real eigenvalues 

Ai < A2 < A3 < • • • ^ +00, 

and there is a unique eigenfunction ipn corresponding to each eigenvalue A„ , which 
has exactly n~l zeros in (0, 1). Moreover, {^n} form an orthonormal base of L|) ^-^ 
under the boundary condition (^(0) = ^{1) = 0. Thus 

Ai = mm- 

V {if, ip) 

where ( , ) denotes inner product. In the case of Theorem 13.21 we also know 
that A2 > by the following argument. Suppose otherwise, A2 < and let ip2 
be its corresponding eigenfunction. Then ip2 has one zero z G (0,1). Let (j) = 
U{y)-U (1/2), then 

(7.1) -0" + g0 = O. 
Also 

(7.2) -¥'2+<3¥'2==A2(^2. 

Multiplying ()7.2p by (p and ()7.ip by (p2, then subtracting, we have 

(7.3) - ip24> + 't'"^2 = A2<^20- 

If z < 1/2, we use the interval [0,z]; otherwise we use the interval [z, 1]. Without 
loss of generality, we assume z < 1/2. Integrating (|7.3p over the interval [0, z], we 
have 

(7.4) - (^^0) |g = A2 / ip2Hy- 



On the interval [0, z], neither (j> nor (p2 changes sign. We note that f2{0) has the 
same sign with ip2 on [0, z], while f2iz) has the opposite sign. Then the right hand 
side and the left hand side of ()7.4|) have opposite signs. A contradiction is reached 
(even in the case (P2{0) — and/or (p2{z) = 0). Thus A2 > 0. □ 

The following lemma is used in the proof of Theorem 16.11 
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Lemma 7.1. Given (/i,/2) G [L"^ {Q.2)Y and J /j^^ /i dydz — 0. Then there 
exists a unique vector field u — {v, w) E (H^ (f^s)) such that 

(7.5) diYu = Vy + Wz = fi, 

(7.6) curlu = Wy -~ Vz ^ f2, 

(7.7) * u-dl^ f u-dl^O, 
and 

(7.8) v{0,z) = v{l,z)^0. 
Moreover, 

(7.9) <C(||/i||^,„ + ||/2||^™) 
for any m > 0. 

Proof. First, we show the existence. We spht u — ui + U2, where ui and 
U2 satisfy 

divui — 0, curl Ml — /2, 

and 

divu2 — fi, curlu2 = 0, 

respectively, with conditions (|7.7p and (|7.8[) . For existence and uniqueness of ui, 
we refer to Theorem 2.2 of [35| Chapter 1]. We construct U2 = V(^, where ip solves 
the Neumann problem 

dtp 

Aip = fi in 512, TT- = on dri2- 
on 

The solvability is ensured by the zero mean assumption on /i and Fredholm al- 
ternative principle. Note that condition (j7.7p is automatic for the gradient flow 

U2- 

Second, we show uniqueness. If there exists two vectors fields u and u' satisfying 
()7.5p - (|7.8p .Then their difference u" = u — u' is an irrotational and divergence-free 
field, tangent to 9i72 with zero circulations on dft2- By Theorem 2.1 of [351 Chapter 
1], u" = and thus u = u' . 

It remains to prove the estimate (|7.9p . For this, we use the following general 
estimate (see (1.26) in |501 Page 318]): Let be an open bounded domain of R^, 
m is an integer > 1. Then for any u G iJ™ (fi) with u • n = on dfl, there exists 
c (m, fi) such that 

(7.10) \\u\\H'r^ < C(||divu||^„_i + ||curlu||^„_i + ||?I||^2). 
We shall show that: for any u e (H^ (fia)) satisfying 1^ and (ITTa . 

(7-11) MIl^q^) - ^0 (l|divw||i2(02) + l|curlM||^2(02)) • 

Then (|7.9p is obvious from (|7.10p and (|7.1ip . We prove (|7.1ip by a contradiction 
argument. Suppose otherwise, for any n > 1, there exists Un G {H^ (^2)) such 
that 

(7-12) NnlL^cn^) ^ (lldivu«|li2(f^^) + |lcurlu„||^2(02)) • 
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we normalize ||m„||^2(s-;2) ^ ^' Then by (I7.10p . ||Mn|| jifi(n2) is uniformly bounded. So 
Un converges to Moo weakly in (il2) and strongly in {^2)- Thus ||'"oo|li2(Q2') — 
1. But from (|7.12p . Uqo is irrotational and divergence- free. Moreover, Uoo also 
satisfies (|7.7p and (|7.8p . So Hoc = 0, a contradiction. This finishes the proof. □ 
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